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The inhomogeneities associated with massless Kalb-Ramond axions can be amplified not only in 
isotropic (four-dimensional) string cosmological models but also in the fully anisotropic case. If the 
background geometry is isotropic, the axions (which are not part of the homogeneous background) 
develop, outside the horizon, growing modes leading, ultimately, to logarithmic energy spectra 
which are "red" in frequency and increase at large distance scales. We show that this conclusion 
can be evaded not only in the case of higher dimensional backgrounds with contracting internal 
qq " dimensions but also in the case of string cosmological scenarios which are completely anisotropic in 

■ four dimensions. In this case the logarithmic energy spectra turn out to be "blue" in frequency 

' and, consequently, decreasing at large distance scales. We elaborate on anisotropic dilaton-driven 

models and we argue that, incidentally, the background models leading to (or flat) logarithmic 
energy spectra for axionic fluctuations are likely to be isotropized by the effect of string tension 
' corrections. 
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In recent years a lot of effort has been devoted to the analysis of rather peculiar phenomenological implications 
■ of string inspired cosmological models [0 . One of the distinctive features of the low-energy string effective action is 
indeed the presence of the dilaton field whose sharp growth amplifies, via gravitational instability, not only the scalar 
and tensor fluctuations of the geometry Q but also the (quantum mechanical) inhomogeneities associated with other 
fields which are not part of the homogeneous background like the gauge fields ||]4j and the universal axion of string 
theory H|), i- e - the (four-dimensional) dual of the Kalb-Ramond antisymmetric tensor field present in the (effective) 
low energy limit of the action. 

The main difference between the amplification of the metric fluctuations and the amplification of other fields with 
no homogeneous background turned out to be, a posteriori, the spectral amplitudes. In the case of metric fluctuations 
the growing modes appearing in the corresponding evolution equations lead typically, to increasing ( "violet" ) spectra, 
whereas in the case of Abelian gauge fields and axions the spectral slope can be much milder. 

These results refer to the case of a homogeneous background geometry whose spatial section naturally decomposes 
into the direct product of two (maximally symmetric) Euclidean sub-manifolds which will be called, for short, external 
and internal. This type of metric describes a situation of dimensional decoupling where the external dimensions expand 
k> , and the internal dimensions contract []. In previous studies concerning the phenomenological implications of string 
inspired scenarios, the expanding dimensions have been always taken to be isotropic (i.e. described by a unique 
scale factor). Similarly, the internal dimensions have been also taken to be isotropic (this situation can happen if we 
compactify the internal dimensions on a six-dimensional torus). This choice is not a limitation. In fact we know that 
the dilaton "vacuum" solutions can be solved for arbitrary (anisotropic) manifolds. 

There are, a priori, no reasons why the external and/or internal sub-manifolds should be considered, separately, to 
be isotropic. In the context of the pre-big-bang scenario [Q the dilaton driven phase starts when the coupling constant 
and the curvature are quite small but this does not necessarily implies the isotropy of the expanding manifold. In 
fact, the study of the occurrence of a dilaton-driven, expanding, regime (in the string frame) can be connected to the 
roblem of the gravitational collapse of a stiff (perfect) fluid in general relativity. Therefore, as it was recently shown 

an anisotropic dilaton-driven phase can be expected also in four dimensions. 
If the four-dimensional metric is not isotropic different phenomenological consequences can be expected both for the 
metric inhomogeneities and for the fluctuations of other fields (like the Kalb-Ramond axions) not contributing to the 
homogeneous background. Concerning the metric perturbations, it is well known that if the four-dimensional geometry 
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is completely anisotropic the scalar, vector and tensor modes of the corresponding fluctuations are all coupled [0. 
More precisely the tensor modes can still be decoupled from the scalar ones, provided the four-dimensional metric 
possess some special symmetry M (like the spherical symmetry, as assumed in Q or the cylindrical symmetry as 
suggested in flOfl ). Similar problems were addressed not in the framework of the pre-big-bang evolution but in the 
context of a post-big-bang (Kasner) behavior of the general relativistic solutions in the vicinity of a cosmological 
singularity ||ll| , |l2| . 

In this paper we will concentrate our attention on the description of mass-less string axions fields in a fully anisotropic 
four-dimensional metric. Dilaton and graviton production in anisotropic string cosmological models can be studied 
using techniques similar to the ones we are going to exploit. There is a simple reason for this choice. In the study 
of the quantum mechanical inhomogeneities amplified in string cosmology the scalar and tensor modes of the metric 
become relevant at large frequencies Q (i.e. small distance scales), whereas the axion and gauge perturbations seem 
to be more important at small frequencies where they can have interesting effects on the cosmic structures. Gauge 
field fluctuations can have interesting implications in the context of the problem of large scale magnetic fields 
and in the context of large scale anisotropies Massless (Kalb-Ramond) axions interact only with couplings of 
gravitational strength and, therefore, their presence is not constrained by present tests of the equivalence principle. 
Moreover the smallness of the coupling to the anomaly makes their conversion to photons (in strong magnetic fields) 
negligible. Recently it was pointed out that Kalb-Ramond axions, amplified from vacuum fluctuations, can repre- 
sent an interesting constraint for pre-big-bang models since they could be invoked for the generation of the CMBR 
anisotropy in string cosmology 

The purpose of our investigation is to understand which are the main consequences of treating axionic perturbations 
in a fully anisotropic four-dimensional metric. Our aim is to show that the evolution of the growing modes of the 
axionic seeds will be significantly influenced by the breaking of the isotropy with quite unexpected consequences. 

Indeed in the case of four dimensional (isotropic) string cosmological metrics we have that the frequency dependence 
of the (logarithmic) energy spectra associated with the axion is typically "red" , i.e. the energy density of axions (in 
critical units) increases at large distances. This occurrence simply indicates that if we wait long enough time these 
fluctuations may become dominant (possibly over-closing the Universe) in the far future. If we go to the case of ten 
dimensional exact solutions of the low energy beta functions (of the type discussed in |^,|^] ) with four expanding and 
six contracting dimensions flat spectra of axionic fluctuations are allowed J5|,^| (provided the axionic fluctuations along 
the internal dimensions are frozen). 

The question we want to address is then the following: if we consider the completely anisotropic case are there any 
indications (already in four dimension) for different spectral behaviors? In this sense our investigation complements 
and extends previous studies on the subject. 

In order to study this problem it is useful to review very briefly the standard results. Consider a four dimensional 
(isotropic) metric of pre-big-bang type |l]]. Then, the Kasner- like nature of the solutions implies that the curvature 
scale and the dilaton coupling arc both growing at the same rate: 



a(r?) = 



V3 + 1 



(1.1) 



where (f> is the dilaton field, and r\ is the conformal time coordinate. In this background the (canonically normalized) 
axionic fluctuation evolve according to 
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C k = 0, C fe = £Vfc, G = ae* 



(1.2) 



( ip is the axion field and the prime denotes derivation with respect to conformal time coordinate; k denotes the 
comoving momentum). Substituting into Eq. (1.2) the background given in Eq. (1.1) we find that the solutions of 
Eq. ( |l.2| ) can be expressed as linear combination of Bessel functions Jl^] with index fi = \/3. It can be checked that 
in the limit when a given mode is outside the horizon (i.e. 77 — >■ 0_ ) the small argument limit of the Bessel functions 
leads to a growing and to a decreasing mode solutions 



Ck(r]) ~ ci(fc) 



V_ 



c 2 (k) 



m 



-V3+i 



(1.3) 



(ci(fc) and C2(k) are two integration constants fixed by the quantum mechanical normalization of the canonical mode 
function). 

If the dilaton driven phase is followed by a radiation dominated phase it is possible to compute the logarithmic 
energy spectrum according to the standard techniques developed in the case of metric and gauge fluctuations 
with the result that the axion energy density will be given, in critical units, as 
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(1.4) 



(u> = k/a is the physical momentum; io\ = k\ja ~ (r/ia) _1 is the maximal amplified frequency; g\ — 0.1 — 0.01 is 
the value of the coupling constant at the moment of the transition from the dilaton-driven phase to the radiation 
dominated phase, and fi 7 (f?) is the critical fraction of energy density stored in radiation at a given time r\). From Eq. 
(1.4) we see that if \i = \/3, as implied by the background (fbl|), the spectral slope in Eq. (1.4) is 3 — 2/i = —0.46, 



namely the logarithmic energy spectrum is decreasing in frequency. Usually these type of spectra are called red simply 
because the largest power is stored at large distance scales. If one wants to use axionic seeds for the explanation 
of large scale (temperature) fluctuations in the Cosmic Microwave Background Radiation (CMBR) red spectra do 
not seem ideal || since one would need either flat or "blue" spectra (i.e. 3 — 2p>0). As noted in |||| by adding 
six (internal) contracting dimensions and by assuming that axion fluctuations do not propagate along the internal 
dimensions, it is possible to achieve \x = 3/2. In this case the scale factors of the external, 3-dimensional (expanding) 
and internal, 6-dimensional (contracting) manifolds will evolve, respectively as a e (rj) ~ 9 — f?) -1 ^ 4 and a;(??) ~ (— rj) 1 ^ 4 
while the dilaton is constant This kind of background could also be dynamically justified not only in terms 

of vacuum solutions but also in terms of a string driven dynamics where the (vacuum) background equations are 
supplemented by (effective) string sources with negative pressure Q . 

Are the conclusions of the present analysis valid in the case of a completely anisotropic four- dimensional metric? 
This is the main point we would like to investigate. 

Thus, we should study which is the evolution of the (super- horizon) axion fluctuations when we drop the assumption 
of isotropy for the four dimensional metric. We will also have to compute precisely the spectral energy density when 
the axionic inhomogeneities, amplified during the anisotropic (dilaton-driven) phase, re-enter during an isotropic 
radiation dominated phase. 

The plan of our paper in then the following. In Section II we will review the basic equations describing the dynamics 
of anisotropic dilaton-driven phases. We will also set up the main evolution equations of the axionic mode functions. 
In Section III we will study the evolution of axionic perturbations in the anisotropic models and we will show, in some 
analytic examples, how blue energy spectra appear. In Section IV we will concentrate our attention on numerical 
solutions. Finally in Section V we will study the influence of string tension corrections on our picture paying special 
attention to the class of models leading to blue axion (logarithmic) energy spectra. Section VI contains some final 
remarks and speculations. 



II. BASIC EQUATIONS 



The low energy string theory effective action in a four-dimensional background is 
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(2.1) 



where H^ va is the antisymmetric tensor field and g^ is a four-dimensional (spatially flat) anisotropic metric 

9fiV = diag[l, -a 2 (t), -b 2 (t), ~c 2 (t)] . (2.2) 

We assume the internal compactification radii to be frozen since we want to analyse, specifically, the breaking of 
isotropy in the external manifold. In order to simplify our discussion we will also restrict our attention to the case 
where b(t) = c(t). This restriction is not essential and all the results discussed in the present paper can be extended to 
the case of b(t) ^ c(t). At the same time we find this simplification useful for the physical intuition. Notice, moreover, 
that if b{t) — c(t) the tensor modes of the geometry propagating along the x axis can be decoupled from the scalar 
modes (this is not the case if b(t) ^ c(t)) ||9j. 

The components of the antisymmetric tensor field JJ_^ va can be re- written in terms o a pseudo-scalar axion as 



(2.3) 



With this choice the action (2.1) becomes 
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(2.4) 
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We want to study the situation where the axion field is not a source of the background but its fluctuations get excited 
by the coupled evolution of the metric and of the dilaton. By varying the action with respect to the dilaton field and 
with respect to the metric we obtain, after linear combinations, the tree-level evolution equations Q for an anisotropic 
string cosmological model 



2cj) + H 2 + 2F 2 = 



H 2 - 2F 2 







H = H(f), F = F(f> 



(2.5) 



where 4> — 4> — H — 2F and the over-dot denotes derivation with respect to the cosmic time coordinate t. The dilaton 
driven (vacuum) solutions of this system can be written as 



o(t) 



b(t) 



<j>(t) = (a + 2(3-1) log 



(2.6) 



with 
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(2.7) 



This last equation implies a Kasner-like relation for the exponents (notice that in the true Kasner case we also would 
have the con dition a + 2/3 — 1 which is not mandatory, in our case, in order to have a solution of the system given 
in Eqs. (|2.5| )). Notice that the solutions given in Eq. ( |2.5| ) have in general two branches connected by scale factor 
duality ]14fl. We will be mainly concerned, in our analysis, with the dilaton-driven branch where both the scale factors 
expand and we will comment, when required, about other (possibly relevant) phenomenological situations. 

For future convenience we can express the actual solutions in conformal time coordinate (related to the cosmic time 
coordinate t by the usual relation a(rj)drj = dt) : 
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By varying the action (2^4) with respect to ip we can obtain the evolution equation of the axionic inhomogeneities 
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which can be also translated in conformal time 



i'" + (0' + 2F)fl>' - V 2 x tp - 
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(2.10) 



where T = (log b)' . Defining the canonical field as C = e^l 2 b ip and going to Fourier space we can further modify Eq. 

( fool ) as 



k 2 +k 2a -- Q - 



C k = 0, G = be? 



(2-11) 



( is the modulus of the longitudinal momentum and kx = yk 2 + k 2 is modulus of the transverse momentum). It 
is clear that this equation reduces to Eq. (L2) if the background is completely isotropic. In fact taking a — ► b and 



recalling that k = \Jk 2 L + k 2 ^, we have that Eq. (1.2) is recovered. At the same time we can clearly see that if a ^ b 



and the background is not isotropic Eq. (1.2) differs qualitatively from Eq. (2.11) which can also be written in a 
slightly different form in terms of k and of the ratio of the transverse and longitudinal momentum r — kx/k^: 
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(2.12) 



In terms of the rapidity y — (1/2) log [(k + /si)/(fc — fcj,)], r = (sinhy) -1 . This form of the evolution equation will 
be useful in Section II for the study of the axi onic ( cano nical) mod e fun ction Ck outside the horizon. Inserting the 
anisotropic dilaton-driven solution of Eqs. (|/7|) and (§!]) into Eq. ( ^7TT| ) we obtain 
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2/x=|2A-l|, A = 



>C k = 0, 



a + 4/3 - 1 
2(1 -a) 



(2.13) 



(2.14) 



Before dealing with the analytic al an d numeri cal s olutions of Eq. ( [2.14 ) we want to present some more qualitative 
considerations. Comparing Eq. ( |l.2| ) to Eq. ( 2.13 ) we can notice some analogies but also crucial differences. Both 
equations can be seen as " Schr oedinger-like" equations but with different "potentials" . More specifically, one can find 
a similarity between Eq. ( p.l3[ ) and the radial part of the Schrocdinger equation in a spherically symmetric potential 
where the term going as r\~ 2 would correspond to the (orbital) angular momentum. The term containing fcy would 
correspond instead to the (central) potential [ER|6|. In our problem there are two relevant cases: either 7 < —2 or 
7 > —2. If 7 < —2, then the "potential" term might become larger, in the limit r\ — ► 0_ (i.e. axion fluctuations 
outside the horizon), than the "angular momentum" term. We can immediately see, amusingly enough, that this 
situation is never realized for the Kasner-like solutions defined by Eq. (^jj). In fact the case 7 < — 2 corresponds, in 
the (a, /3) plane to /3 > 1 which is not realized since in the (anisotropic) vacuum solutions a and /3 have to lie on the 
ellipse a 2 + 2/3 2 = 1, and, therefore we will always have that |/3| < 1/V2 ~ 0.707.... Thus, we could always expect the 
r\~ 2 term to be, at some point, dominant for rj — > 0_. 




( [2,g| ) and the behavior of the axionic mode 
All the a and (5 on the ellipse give rise to 
2D with c = b. The arc bounded by the two 



FIG. 1. We illustrate the interplay between the solutions given in Eqs. (2 
function outside the horizon. The ellipse is the geometrical version of Eq. (2 
solutions of the low energy beta functions in the anisotropic metric of Eq. (\ 
full (thick) lines defines the class of anisotropic models leading to /1 < 3/2. The shading is only meant to guide the eye. The 
intersections of the full lines with the ellipse defi ne th e models for which n — 3/2. The dashed curve correspond to /3 = a. For 
j3 > a (a bove the dashed line) we have, in Eq. (2.13) that 7 < 0. The dot-dashed line corresponds to the case where A = 1/2 
(see Eq. (2.14)), implying u — 0. 



In order to visu alize our qualitative arguments, let us look at Fig. |l| wh ere the ellipse is just the plot (in the (a, 0) 
plane) of Eq. (2.7) defining the Kasner-like solution of Eqs. (2.5) and ( |2.8| ). The dashed line corresponds to j3 = a. If 
(3 > a (above the dashed line) 7 < and the "potential" term goes to zero for rj — > —00. Since rf 2 becomes dominant 
for rj — > 0_, the solutions of Eq. ( |2 . 13[ ) should always go to Bessel functions with index p,. 



In the anisotropic case the index [i can be different from the Bessel index appearing in Eq. (1.2) for the isotropic 
case. Given that the (logarithmic) energy spectrum is determined by the dominant solution of Eq. (2.13) in the 
limit rj — > 0_, we can expect for /Lt<53/2 a flat or bl ue log arithmic energy spectrum. In the (a, 0) plane the region 
corresponding to /i<3/2 can be obtained from Eq. ( 2.14 ) with some elementary algebra. The result is reported in 
Fig. |l|. The arc of the ellipse bounded by the two (full) thick lines correspond to (anisotropic) vacuum solutions 
leading to growing modes compatible with "blue" spectra. The intersections between the two full lines and the ellipse 
correspond to the case n = 3/2 (flat spectrum). Of course not all the models within the full lines are realistic. 

For examples we could have, in principle, growing spectra with a > and (3 > (first quadrant). These models are 
clearly not realistic since they would describe a dilaton-driven contraction more than a dilaton-driven expansion. More 
realistic are the models in the second quadrant where one of the two scale factors expands and the other contracts. 

Finally in the third quadrant we have the most physical situation where the two scale factors both expand, though 
at different rates. It is amusing to notice that the intersection between the full line and the ellipse picks up a dilaton- 
driven solution with pi = 3/2 and expanding scale factors. The lower full line has e quat ion 4/3 = (a — 1) and the 
intersection we are interested occurs for a = —7/9 and j3 — —4/9. Going to Eq. (2.8) we see that, in this case 
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4>{rf) ~ —3/2 log [ — 77/ 771] . In this case, on the basis of the asymptotic behavior of Eq. (2.13) in the small \rj\ limit, the 
solution outside of the horizon will be given by 



Ck(v) — ci(fc) 



c 2 (k) 



m 



(2.15) 



which coincides with the solution one would obtain in the ten-dimensional case where the external manifold has scale 
factor a e (r]) ~ (— ^) -1 / 4 and the internal one has scale factor ~ (— 77) 1 / 4 . 

Another interesting case is represented by the model a = —1 and [3 — 0. In this case 4>(rj) ~ — log [—77/771]. This 
solution might not seem fully realistic and indeed we will use it as a toy model. At the same time one can argue that 
this type of solutions might be the result, in the pre-big-bang scenario, of spherically symmetric initial conditions 
which are asymptotically trivial in the past [Q . It is quite amusing to work out the Einstein frame description of this 
solution. The Einstein frame scale factors are related to the ones of the String frame by trivial conformal rescalings 
involving the dilaton, namely 



b = e 2 b. 



(2.16) 



Since the conformal time coordinate does not change passing from one frame to the other (i.e. drj = df) we can obtain 
easily that, in the Einstein frame, our solution looks like a ~ const, and b ~ (since a(rj) is constant t and 77 

coincide up to a constant factor). 

It is interesting to point out that this solution was Kasner-like in the String frame (i.e. a 2 + 2/3 2 = 1) and it becomes 
true Kasner in the Einstein frame. This solution is well known in the theory of the gravitational collapse jL7| and it 
might be thought as the limit (in the vicinity of the singularity) of a spherically symmetric space filled with a stiff 
(perfect) fluid. It was indeed recently stressed that the Einstein frame description dilaton-driven cosmologies in four 
dimensions is closely connected with the problem of the gravitational collapse of a stiff fluid and that, consequently, 
string cosmological models of pre-big-bang t ype a re as fine-tuned as the gravitational collapse Q] . 

If a = — 1 and j3 = we have that in Eq. (2.13) 7 = — 1. It is amusing to notice that in this case the Eq. (2.13) is 



analytically solvable and it is (formally) equivalent to the Schroedinger equation in a Coulomb potential with complex 
electromagnetic coupling. We w ill us e this analogy in the next section. For the moment we can guess that outside 
the horizon the solution to Eq. (2.13) will be well approximated by 



C k = ci(fc) 
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(2.17) 



Again this solution was simply obtained by solving Eq. ( 2.13| ) in the small 77 limit where the "Coulomb" contribution 
can be neglected. Since this is nothing but the small argument limit of Hankcl functions with ji = 1, we again argue 
that the exact solution should give us a blue (logarithmic) energy spectrum of the type of the one given in Eq. ( |l.2| ) 
with spectral index 3 — 2/j = 1. 



III. BLUE SPECTRA 



In some cases, analytical solutions of Eq. ( |2.13| ) c&n be found. We start to discuss blu e spe ctra and, more precisely 
we analyze the model mentioned in the previous Section with a = —1. In this case Eq. (2.13) becomes, using k and r 
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This equat ion looks indeed very similar to the (radial) Schroedinger equation for the Coulomb problem. In order to 
solve Eq. (3J) we have to match the standard notations appearing in the theory of special functions [(TBI . Defining a 
new (complex) time as 



2ikr) 



we can re- write Eq. (3.1) as 
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with £ = [i/2)er. It is helpful to define e = kxTji = krjir/y/l + r 2 purely for algebraic reasons. Eq. ( |3.3| ) is one of 
the well known forms of the Whittaker's equation whose solutions are the corresponding Whittaker's functions, 
related to confluent hypergeometric functions. We are interested in solutions which represent incoming axionic waves 
for rj — > — oo and then our solution will have the form 



W h {r)=e-*T*U\ 



(3.4) 



where U[m, n,r] is the (Kummer) confluent h ype rgeometric function with parameters m and n. Therefore the (prop- 
erly normalized) form of the solution of Eq. (3.1) is 
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(3.5) 



We can immediately notice, as expected, that the sma ll arg ument limit of the previous solution (i.e. kr\ <C 1) exactly 
reproduces the asymptotic evolution o btai ned in Eq. (2.17). 

Having derived the solution of Eq. (3.1) we can also compute the energy spectrum of the axions amplified by the 
transition from an anisotropic dilaton-driven solution to a (completely isotropic) radiation dominated epoch occurring 
at 77 = —T)i. For r\ > —rji the evolution equation of the axion fluctuations will be given by 



C'l + k 2 C k = 0, 



(3.6) 



because the dilaton coupling freezes, after —771, to a constant value. A similar calculation was performed in Q under 
the assumption that the isotropisation of the geometry occurred prior to the onset of the dilaton-driven phase (indeed 
described, in || only by one scale factor). 

The frequency mixing coefficient c_(/c), determining the spectral numb er o f produced axions, can be explicitly 
computed by matching, at 7/ = —rji the general "outgoing" solution of Eq. (3.6) 



Cfefa) = 



c+e 



-ifc(?j+r/i 
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V > -vi, 



(3.7) 



to the "incoming" solution ( |3.5| ) of Eq. ( |3.l| ) which represents positive frequency modes in the rj — > —00 limit (notice, 
however, that the incoming solution does not define, asymptotically a Bunch-Davies vacuum because of the presence 
of the term "Coulomb" term in Eq. 



(3.1 



which is dominant in the far past). 
Following this procedure we determined the mixing coefficients reported in Appendix A. The important point to 
mention is that for krji > 1 the procedure le adin g to the mixing coefficients is no longer valid but the mixing coefficient 
can be anyway estimated by using in Eq. (3.1) a smooth "potential" interpolating between 77 < — 771 and 7/ > —7/1. 
For krji > 1 the mixing coefficients are indeed exponentially suppressed [ fl8| and can be ignored for the purpose of 
this paper. Therefore, taking the limit (for krji < 1) of the results reported in Appendix A we have 



[S(r)) 



-1/2 



-3/2 



47T 



(3.8) 



where S(r) = 1/ \Jr 2 + 1 and x\ = kr\\. 
The spectral energy density 



dp^p 



n . „ (3.9) 

CI lOg LO 7T Z 

is the variable frequently adopted to characterize the distribution of the produced particles dHII > f° r our case rea ds 

gf0 7 (7?) ( u; 



77) 



(3.10) 



where p c is the critical energy density and (a;, 77) is the (critical) fraction of produced (mass-less) Kalb-Ramond 
axions. Notice that n{u) = |c_(cj)| 2 can be interpreted, in a fully quantum mechanical treatment of the pro cesse s 
of cosmological particle production P,[l9|, as the mean number of axions with physical momentum uj. In Eq. ( 3.K ) 
gi = e"^'' 1 '/ 2 ~ 0.1 — 0.01 is the final value of the string coupling parameter at the time of the transition to the 
radiation dominated epoch. In terms of gi the string and Planck mass are related as M s = g\Mp. 
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The s pectr al slope appearing in Eq. (3.10) is exactly the one we guessed in Section II using the approximate solution 
of Eq. ( 2.13 ) in the small 77 limit for the axion mode function. The result is that the spectral energy density is only 
slightly tilted towards large frequencies. We remind, incidentalfythat in string cosmological models blue spectra are 
also possible for gauge fields 0], but not for metric fluctuations g. In the case of gravitons, for example, the typical 
slope for the spectral energy density is close to 3 ("violet" energy spectra). 

Thus also in the fully anisotropic case the spectral distribution of the amplified axions is determined by the solution 
of the mode equation in the small 77 limit. Furthermore, on the b asis of the analytical treatment performed in the 
present Section we can argue that as long as 7 > —2 in Eq. ( 2.13 ), the spectral slop e of t he spectral energy density 
p u is completely determined by the "angular momentum" term appearing in Eq. (2.13). On the other hand the 
spectral amplitude does depend upon the details of the particular anisotropic model and upon the relative weight of 
the longitudinal and transverse momentum (i.e. r = kx/k^) during the anisotropic phase. In particular if r — > (i.e. 
kx <C fcj,) S(r) — ► 1. In the large r limit (i.e. kx 3> fcz,) S(r) ~ r. From a physical point of view, r ~ 1 implies a fairly 
isotropic (initial) momentum distribution and it makes sense only when a(r[) — > b(rj). Conversely, r very different 
from 1 implies either the dominance of transverse momenta or the dominance of the longitudinal momenta. We see 
that r cannot be arbitrarily large since we have always to demand fl(u;,r)) < 1. Imposing the critical energy bound 
on the axion fluctuations we get, indeed, that r < I6ir 4 /gf (i.e. r < 10 5 if g x = 0.1). 



One can understand, more physically, the r and krji dependence by looking more closely at Eq. (3d). Let us write 
the "Coulomb" potential and the "angular momentum" in a rescaled form, namely 
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where £ = [— 77/771] and, as usual, x\ — krj\. The plot of the function (3.11) for (different values of the parameters) is 
reported in Fig. 
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FIG. 2. We plot different regimes of Eq. (3.11). With the thick (full) line we have the case x\ = 1 and r = 1. With the 
dashed line we illustrate the case X\ — 0.1, r = 1, whereas with the dot-dashed line we illustrate the case x\ — 1 and r = 10 3 . 
Notice that the thin (full) line is for the case x\ = 0.001, r — 1. 



Since we want to consider modes which are amplified (and not exponentially suppressed) we have to focus on the 
case xi<Jl. From Fig. ^| we learn that as soon as x\ gets smaller and smaller (i.e. frequencies smaller than the 
maximal) the "Coulomb" term in Eq. ( |3.11 ) gets even more negligible if compared with the term going as £ -2 . 
Moreover, for large r and high frequencies (i.e. x\ ~ 1) the potential barrier gets higher. Thus, if we are going to 
explore a frequency range k -C 77^ 1 (i.e. uj <C u>x) our approximation improves. 



IV. FLAT SPECTRA 

Having discussed the case of a specific blue spectrum we want now to treat the case of a flat spectrum in anisotropic 
dilaton-driven models. We argued in Section II that in order to obtain a flat spectrum we have to require /1 = 3/2, 
which does correspond either to a = 1 and f3 = or to a — —7/9 and (3 = —4/9. The case a = 1 (3 — is clearly not 
realistic : it does correspond to a contracting Universe (in the string frame) . The interesting case is indeed the second 
where the dilaton growth triggers a superinflationary phase where the different dimensions expand with different rates: 
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For this background the explicit for of Eq. (2.13) is 
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which becomes, going to dimensionless (conformal) time, 
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I n E q. (4.4) the analogous of the Coulomb term goes as g~ 3 / 8 . Therefore, we expect the (numerical) solutions of Eq. 
(4.3) to converge more slowly than in the case of Eq. ( 2.17 ) to the solution one would guess taking into account only 
the leading contribution for small £. This is illustrated in Fig. || where the function is reported for different 
values of the parameters. 
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FIG. 3. We illustrate different cases of Eq. (4.4). With the thick (full) line we have the case Xi 
dashed line we plot the case xi = 0.1, r = 10. 



1 and r = 1. With the 



Eq. (4.3) can be solved numerically. In Fig. (||) the numerical solutions of Eq. (4.3) are reported for different 
values of the primordial rapidity. For comparison (d ashe d lines) we also report, for the same values of the parameters, 
the behavior of the approximate solutions of Eq. (4.3) obtained by neglecting the term £~ 3 / 8 . For £!<(). 001 the 
approximate solutions can be definitely trusted in the small £ limit. 




FIG. 4. We report the numerical sol utio ns of Eq. (4.3) for decreasing values of xi. We see also report the approximate 
solutions (dashed curves) given in Eq. (4.5). 







The approximate solutions plotted with the dashed line in Fig. (0) are 
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By directly matching Eq. ( |4.q ) with Eq. (3.7) we obtain that the relevant mixing coefficient is 
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corresponding to 
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Notice that fC{r) is a (frequency independent) factor which can be precisely determined numerically and which becomes 
important in the limit x\ ~ 1. If xi<0.001, IC(r) ~ 1. 

It would be desirable to find approximate (but analytical) expressions describing the evolution of the mode fun ctio n 
outside of the horizon in the general case of anisotropic backgrounds of the type of the one introduced in Eq. ( |2.S| ) . 
The evolution equation (2.11) for the axionic mode function can be written as 
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The effective potential appearing in Eq. (4.8) generalizes the previous results of Eqs. ( 3.11 ) and (4.4). Notice that 
the situation of a fc-dependent effective potential is not so uncommon since it arises in the case of perturbations of 
kaluza-Klein backgrounds with fluctuating internal dimensions (|] and in the case gravitational waves propagating 
in a background whose action contains higher derivatives Notice also that in the isotropic limit (i.e. b — > a) 

Eq. (4.8) goes to Eq. (1.2) since Q — > ae"^ 2 and the term proportional to the anisotropy (i.e. (b 2 — a 2 )/b 2 ) exactly 
vanishes. We want to solve Eq. (4.8) in the limit — ^ 0_ and show that the correction to the leading order solution 
(obtained by dropping the b 2 — a 2 )/b 2 term) is indeed small. The equation we have to solve is essentially 
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By writing the solution to Eq 
where Ck is the leading order solution 



C'l ~ V k (r,)C k . 
Ck{ri)~Ck{ri) + ek{ri), 

C k ~ AG( V ) + BG( V ) 



drj' 



we can find that e-kijl) is simply given by 

efefa) 



1 



G(V) / G(t)C(t) 



b 2 ~a< 
b 2 



drf_ 
G(v") 2 



(4.9) 
(4.10) 

(4.11) 
(4.12) 



(A and B are arbitrary constants). Notice, again, t hat t he correction vanishes exactly for b —> a. If we insert in the 
previous equations the background solutions of Eq. (2.8) we get 
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In each of the curly brackets the first term is the leading solution whereas the term proportional to (km) 2 is the 

correction. Notice that the correction is subleading in the limit 77 — > Take, for instance, the case investigated in 

this section, namely A = — 1 and 7 = —3/8. We can easily find, from Eq. (4.13) 
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Notice that the leading super-horizon contribution to the evolution of the axionic mode function is given, as expected, 
by the term going as n" 1 appearing in the first curly bracket. In both the curly bracket we have the term (fcr/i) 2 
which is always smaller than one since 77^ 1 is the maximal amplified frequency. Moreover, in both the curly brackets 
all the terms multiplied by (fcryi) 2 are vanishingly small outside the h orizon . We could repeat the same calculation 



by inserting different (expanding) anisotropic backgrounds into Eq. ( 4.15 ). We find that for the class of models 
defined on the arc of the "vacuum" ellipse depicted in Fig. |l| all the corrections are inmdeed small and the leading 
contribution is given, as foreseen, by the growing mode. Notice that we already knew that from the considerations 
of Section I and II. We the present approach we are also able to compute, analytically, the corrections to the leading 
solution. 

V. ANISOTROPY AND STRING TENSION CORRECTION 

The main assumption behind the calculations we presented has been the presence of an anisotropic string cos- 
mological phase. In previous calculations the axion spectra were actually computed assuming that the (four 
dimensional) background describing the dilaton-driven dynamics was completely isotropic. Therefore, in |5||| it was 
also implicitly assumed that the isotropization of the four (expanding) dimensions occurred prior to the onset of the 
dilaton-driven evolution. As we stressed in Section II, to assume that the dilaton driven phase was not completely 
isotropic might be natural from the point of view of the "asymptotic past triviality" of the pre-big-bang scenario . 
More specifically, the comparative study of (spherically symmetric) gravitational collapse and pre-big-bang scenario 
seem to imply the occurrence of anisotropic dilaton driven phases. As we stressed, the example given in Section III can 
be interpreted, from the point of view of the background evolution in the Einstein frame, as the gravitational collapse 
of a sphere filled with a perfect (stiff) fluid ]17fl leading, ultimately, to an anisotropic four-dimensional dilaton-driven 
phase in the string frame picture. 

The fact that anisotropic initial conditions are generally allowed by the solutions of the tree-level action might not 
be, in principle, sufficient in order to justify our assumptions. In fact what we need for our considerations is a metric 
which is anisotropic not only locally (i.e. for a short amount of time in the dilaton-driven phase) but globally (i.e. 
for all the duration of the dilaton-driven phase). What might happen is that the addition of higher string tension 
corrections to the tree-level action forbid the presence of a (long) anisotropic phase. If the anisotropic phase is not 
sufficiently long the axion fluctuations will not have enough time to fully develop their growing modes. 

A fair measure of the anisotropy of a four-dimensional metric can be given, in our particular case as 

W V H + 2F ' W 3 V ; 

(this type of measure for the anisotropy of a four-dimensional geometry was introduced long ago by Zeldovich |po| ) . 
Notice that, in the isotropic limit (i.e. H — ► F), V(t) coincides with the Hubble factor and A(t) — > 0. Another 
possible measure of the degree of anisotropy of a (homogeneous) background geometry can be also given as the ratio 
between the Weyl and Riemann invariants 

B l t ) = (5.2) 

which is supposed to vanish in the isotropic limit (see Appendix B for the explicit expression of B{t) in our background 
geometry) . 

Consider now a dilaton driven-phase evolving according to the solutions discussed in the previous section, i.e. 
H ~ — 7/(9t) and F ~ — 4/(9i). Then A(t) ~ 3/5. The question is then the following: if we introduce quadratic 
corrections to the tree-level action will A(t) be roughly constant for the whole duration of the dilaton driven phase? 
In other words: can we show that A is roughly frozen to the value given by the dilaton-driven "vacuum" solutions 
prior to the onset of the stringy phase where the curvature is stabilized? If this is not the case, namely if A(t) is not 
constant during the dilaton-driven phase, then the anisotropy in the background will only be local. 
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Let us start from the string effective action supplemented by the first string tension correction in four (anisotropic) 
dimensions. We follow the notation of |^2j where this problem was studied in the string frame picture (see also 0] 
for the discussion of similar problems in the Einstein frame picture): 
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where R 2 GB is simply the Gauss-Bonnet invariant expressed in terms of the Riemann, Ricci and scalar curvature 
invariants 
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and w is a numerical constant of order 1 which can be precisely computed depending upon the specific theory we 
deal with (for instance w = —1/8 for heterotic strings). Notice that in jl9| this problem was discussed in greater 
generality. Here we will recall and extend the topics which are directly relevant in our investigation. 

It is convenient, for the problem at hand, not to synchronize immediately the clocks. Therefore the metric will be 



9liv = diag[7V(t) 2 , -a(t) 2 , -6(t) 2 , ~b{tf 



(5.5) 



Notice th at b y fixing the lapse function to 1 corresponds to the synchronous coordinate frame. Inserting Eq. (5.5) 
into Eq. (5.3) we get, after integrating by parts the derivatives of the lapse function, 
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where we used the notation a(t) = e a ^\ b(t) 
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We can take the variation of this action with respect to N, a, 



(3 and <j> and we will get three dynamical equations supplemented by a constraint connecting the first derivatives of 
the fields. These equations are reported in Appendix B. Here we give the final expression of the system of nonlinear 
differential equations obeyed by x = a, y = [3, z = <fi: 
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(5.7) 

(5.8) 
(5.9) 



It is known from previous studies ]19{ that this system of nonlinear differential equations a dmit s an isotropic fixed 
point which can be obtained by setting x — y — constant in the constraint given in Eq. (B.l) which will give a 
consistency condition which leading to x — y — 0.616... and z — 1.414... . It is also known that this (isotropic) 
fixed point can attract (mildly) anisotropic dilaton-driven cosmologies Jl9[ . The question we have to address, in our 
context, is threefold. First of all we have to ask if the anisotropic models we are interested in from the previous 
sections are attracted towards this isotropic fixed point. If this is the case we have to investigate how the anisotropic 
phase is analytically connected with the isotropic phase with constant (and regular) curvature. Finally we have to 
understand quantitatively how long is the duration of the anisotropic phase before the isotropization (triggered by 
the string tension corrections) takes place. 

In our present notation we have, from Eq. (pTTl) that 



A(t) = 



3[*(*) - y(t)} 

x{t) + 2y(t) ' 



(5.10) 



In order to answer the previous questions we will integrate numerically the system given in Eqs. (5.7)-( p.9[ ) imposing, 
at the initial integration time, A(to) — 3/5 which is equivalent to impose initial conditions will along the lines 

4> = (9/7)H and <\> = (9/4)F in the ((f), H) and (<f>,F) planes. Our results are reported in Fig ||. We see, in this 
particular case, that the model, initially anisotropic, is indeed attracted towards its isotropic fixed point by virtue of 
the string tension corrections. 
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FIG. 5. In the two upper plots we report the behavior of x, y and z. At the top left the evolution of z (thick full line) is 
compared with the evolution of x (dot-dashed thick line) . We see that they both freeze to a constant value given by the isotropic 
fixed point. Also y reaches, for large t its isotropic limit but, during the dilaton-driven phase the background is anisotropic. 
In the top right plot we show indeed the evolution of x (thin dot-dashed line) and y (thin full line) prior to the isotropization 
occurring at the onset of the stringy phase. This point is also stressed in the two lower plots. In the bottom left plot we report 
the numerical evaluation of the asymmetry given in Eq. (^l|) . We can see that it goes to zero for large positive times but it is 
roughly constant during the dilaton-driven phase. In the bottom right plo t we i llust rate the same qualitative evolution of the 
ratio between the Weyl and the Riemann invariants reported in Eqs. (p.2|) and (B.7). 



The curvature invariants (inc luding the Weyl invariant) are all regular. Moreover the Weyl invariant and the 
asymmetry defined in Eq. (5.1) both vanish for large times signaling that the isotropy is indeed recovered. The 
anisotropic and the isotropic phase are then connected and the anisotropic regime lasts as long as B{t) and A(t) are 
almost frozen to their constant values (fixed by the initial conditions). This evolution is typical, in our examples, for 
negative times. In this regime the solutions are the ones discussed in Section II and, therefore the typical behavior of 
the axionic growing modes is different from the one dictate d by a com pletely isotropic dilaton driven dynamics and 
leading to red logarithmic energy spectra. By solving Eqs. (5.7)— (5.9) together with Eq. ( 2.1C ), we checked that, as 
far as the dilaton driven phase is concerned, the axionic growing modes are indeed the ones we computed by using 
the tree-level solutions discussed in Section II. We will report this study elsewhere. 

One can ask, at this point, how general is this behavior. Our answer is that it holds for various models which 
we would call mildly anisotropic, from the point of view of their initial (primordial) asymmetry. The results we just 
presented in a particular case do not hold in the case of strongly anisotropic initial conditions where, for example, 
one of the two scale factors contracts (in the string frame) instead of expanding like the ones appearing in the second 
quadrant of the plot reported in Fig. [j], like, for instance the case where a = —1/3 and (3 = 2/3. Since for the 
(anisotropic) dilaton-driven solutions discussed in Section II we have that the asymmetry is 



A(t ) 



>-/?) 
' a + 2(3 



(5.11) 



we have for a = —1/3 and f3 = 2/3 that A(to) = —3. With these initial conditions the system of Eqs. ( 5.7 )— (p79|) 
evolves towards a singularity, the isotropic fixed point is not reached and the anisotropic phase lasts forever. Moreover, 
we can see that the asymmetry grows (im modulus) becoming more and more negative. In this case the string tension 
corrections will make the geometry more and more anisotropic. 

Therefore we studied the acce lera te d br anch of the dilaton driven solutions (i.e. a < and (3 < 0) and we fixed the 
initial conditions of the system ( 5.7)— (5.9) by requiring that 



A(t ) 



' -v/l-2/3 2 + 2/3 



(5.12) 
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(where we simply inserted a = —y/l — 2(3 2 into Eq. ( 5.11 ) in order to force the initial conditions to lie inside the 
space of the solutions of the tree level action [see also Fig. m). We then moved (3 in the range — 0.333(3)</3 < where, 
according to Fig. [l], blue axionic spectra are expected on the basis of the tree-level dilaton driven solutions. Our 
results are reported in Fig. ^ in terms of the two previously introduced estimators of the anisotropy A(t) and B(t). 
We plot four cases namely (3 = —0.333(3), —0.222(2), —0.11(1), —0.55(5) and we find the same behavior discussed in 
Fig. H though with different initial asymmetry. Notice that the model with (3 — —0.555(5) does not give blue spectra 
but, nonetheless, it is driven towards isotropy for large positive times. 
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FIG. 6. Behavior of the asymmetries A(t) and B(t) for different initial conditions within the region of the (a, j3) plane leading 
to blue axion spectra. 



We recall that the motivation of this last section was to understand if string tension corrections can prevent the 
existence of anisotropic phases. In this sense, our conclusion is that string tension corrections do not prevent the 
existence of long anisotropic dilaton-driven epochs. Moreover, we showed that the models leading to blue axion 
spectra are likely to be attracted towards isotropic fixed points as a result of the string tension corrections. 



VI. CONCLUDING REMARKS AND SPECULATIONS 

In this paper we relaxed the isotropy assumption in the study of the amplification of the axionic fluctuations in 
string cosmology. In the isotropic case the four-dimensional spectra are typically red. We then focused our attention 
on the fully anisotropic (four-dimensional) dilaton driven solutions where the growing modes of the axion fluctuations 
are completely different from the one occurring in the completely isotropic case with the same number of dimensions. 
Assuming (as in the isotropic case) a smooth transition to a radiation dominated phase, the obtained axion spectra 
have an amplitude which depends upon the primordial anisotropy of the space-time encoded in the asymmetry between 
longitudinal and transverse momenta which can be expressed in terms of a primordial rapidity. The spectral slope is 
determined by the growing mode solution and can be cither flat or blue for anisotropic models whose scale factors are 
all expanding though at a different rate. Our conclusions are summarized in Fig. m 
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FIG. 7. Axionic logarithmic energy spectra for different dilaton-driven models with Kasner-like exponents (a, (3) in the string 
frame. 

As we showed the situation is qualitatively different if compared with the isotropic dilaton-driven evolution where 
only red spectra are possible. We focused our attention on blue spectra because they are produced by models with 
expanding four-volume and growing dilaton coupling. Red and violet spectra seem less natural from the anisotropic 
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point of view since they require that at least one of the scale factors contracts (in the string frame) . We compared our 
phcnomenological point of view with anisotropic models occurring as a result of spherically symmetric gravitational 
collapse and we showed that some of those models lead to blue axionic energy spectra. This conclusion extends and 
complements previous works on the subject ^,[|. 

Motivated by our phenomenological considerations we investigated the influence of string tension correction on the 
evolution of anisotropic dilaton-driven epochs. We found that string tension corrections never prevent the existence 
of anisotropic phases. Moreover if the anisotropy is mild (like in the region of Fig. |^ compatible with blue spectra) 
all the pre-big bang phase is anisotropic but it is attracted for large positive times towards an isotropic fixed point. 

What we presented in this paper seem to suggest that it might be plausible to consider more seriously the possible 
effect of primordial anisotropics on the axionic spectra as it was recently done in the purely isotropic case ||. It 
is not impossible to imagine that the amplitude of our spectra (and then the primordial ratio between longitudinal 
and transverse momenta) can be constrained by the temperature fluctuations in the Cosmic Microwave Background 
radiation. More generally, we find interesting that our considerations seem to point towards a more detailed analysis 
of anisotropic string cosmological models. 
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APPENDIX A: MIXING COEFFICIENTS 



In this appendix we report the mixing coefficients employed in the calculation of the energy spectrum of the amplified 
axionic perturbations discussed in Section III. Defining the function 

S(r) = -=J=, r= k -L^ (A.l) 

we have that the mixing coefficients are 

M*) = ™^e-^)) { (-3 - Ai X + || - ^(3 + - ler, 3, ||] 

2a; -(3z + 6r)C/[|~^r,4,2z^(r)]} 



S(r) K ' L 2 2 

c+(fc) = ffi^^ e -teP-^)] {( _3 + 4 zx + || + r 2 (-3 + 4*a:))Z7[| - ler, 4, 2isS(r)] 

2S -(3z + 6r)C/[^~^r,4,2z^(r)]} (A.2) 



S{ry ' L 2 2 



where x = krj. 



APPENDIX B: EQUATIONS OF MOTION WITH STRING TENSION CORRECTIONS IN THE FULLY 

ANISOTROPIC CASE 



Taking the variation of the actin given in Eq. (5.6) with respect to the lapse function N(t) and imposing, afterwards, 
the cosmic time gauge we get the constraint 



4? + 2/3 2 + 4a(3 - 2acj> - AP<\> + 



3 ■ 



4 - 6a60 2 
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(B.l) 



where we took string units X s = 1 and w = 1. By varying the action with respect to a, (3 and <f> we get, for N(t) = 1, 



4,9-20 + (4/3 -2(f>) (d + 20- 



(d + 2(3 - 0)0/3 2 + 0/3 2 + 2/3/30 



+ L(t) = 0, 



4(/3 + a - <f>) + 4(d + /3 - 0)(d + 2/3 - 0) - 
2(0 - a - 2/3) + 2(0 - a - 2(3){a + 2/3-0) 
- 20(a0 + 2a.fi) + {2a(3 2 - 3 )(ci + 2/3 - 0) - 30 2 

where we defined 



/3a<j> + $a(j> + $a<j> + /3acb(a + 2(3 - 



2L{t) = 0, 



L(t) = 0, 



Defining now 



L(t) = -0 2 - 2(3 2 - Aa(3 + 2a0 + 4/30 + -(80d/3 2 - 4 ). 



x(t) = a, y(t) = 0, z(t) = 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



we get, after trivial algebra, the equations reported in Section VI. 

In Section VI we also used an alternative measure of the anisotropy of a homogeneous space-time, namely the ratio 
between the Wcyl and Ricmann invariants. In the specific case of our metric this quantity reads: 



B(t) = 



C, 



[xy — x 2 — x + yY 



R^ a f3R^ afs 3 [x A + 2x 2 y 2 + 3y 4 + 2x 2 x + x 2 + Ay 2 ij + 2y 2 } 
We can clearly see from the above expression that in the isotropic limit (i.e. x — y) Bit) — > 0, as expected. 



(B.7) 
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